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Abstract. We present a general formalism for deriving bounds on the shape parameters of the weak and 
electromagnetic form factors using as input correlators calculated from perturbative QCD, and exploiting 
analyticity and unitarity. The values resulting from the symmetries of QCD at low energies or from lattice 
calculations at special points inside the analyticity domain can be included in an exact way. We write down 
the general solution of the corresponding Meiman problem for an arbitrary number of interior constraints 
and the integral equations that allow one to include the phase of the form factor along a part of the 
unitarity cut. A formalism that includes the phase and some information on the modulus along a part of 
the cut is also given. For illustration we present constraints on the slope and curvature of the K13 scalar 
form factor and discuss our findings in some detail. The techniques are useful for checking the consistency 
of various inputs and for controlling the parameterizations of the form factors entering precision predictions 
in flavor physics. 



1 Introduction 

Form factors are of central importance in strong inter- 
action dynamics, providing information on the nature of 
the strong force and confinement. Phenomenologically, the 
weak form factors are of crucial importance for the deter- 
mination of standard model parameters such as the ele- 
ments of the Cabibbo-Kobayashi-Maskawa matrix. 

Bounds on the form factors when a suitable integral of 
the modulus squared along the unitarity cut is known from 
an independent source were considered in the 1970's in the 
context of the hadronic contribution to muon anomaly 
and the kaon semileptonic decays [Il[a|3lll[5l[6j[7l[8l|9l[l0l 
[TTj (for a topical review of the results at that time, see 
[12]). Through complex analysis, this condition leads to 
constraints on the values at interior points or on the ex- 
pansion parameters around t ~ 0, such as the slope and 
curvature, and belong to a class of problems referred to 
as the Meiman problem [13] . Mathematically, the problem 
belongs to the standard analytic interpolation theory for 
functions in the Hardy class [HlHtlllSlliej . The inte- 
gral condition was provided either from an observable (like 
muon's g — 2 in the case of the electromagnetic form factor 
of the pion) , or from the dispersion relation satisfied by a 
suitable correlator, whose positive spectral function has, 
by unitarity, a lower bound involving the modulus squared 
of the relevant form factor. Therefore, the constraints de- 
rived in this framework are often referred to as "unitarity 
bounds" . 

An important step forward was achieved in jl7j , where 
it was noted that the correlator of interest for the K13 form 



factors can be evaluated reliably in the deep Euclidean 
region by perturbative QCD. The modern approach clar- 
ified also the issue of the number of subtractions required 
in the dispersion relation for the correlator, which was not 
always treated correctly in the prior studies. 

The first applications in the modern approach con- 
cerned mainly the form factors relevant for the B ^ D and 
B — )■ D* semileptonic decay, or the so-called Isgur-Wise 
function, where heavy quark symmetry provided strong 
additional constraints at interior points [TSllTOll^ II 2 1 ll^ 
[23l[24l[25] . More recent applications revisited the electro- 
magnetic form factor of the pion [26 . 27"281l29j , the strangeness 
changing Ktt form factors {501151 . 32. ,33^ . and the Btt vec- 
tor form factor p4ll35j . The results confirm that the ap- 
proach represents a useful tool in the study of the form 
factors, complementary and free of additional assumptions 
inherent in standard dispersion relations. 

The purpose of the present paper is to present in a sys- 
tematic way the technique of unitarity bounds, including 
its most recent developments and offering explicit formu- 
las that can be applied in future studies. There arc several 
recent developments which increase the interest in these 
techniques and justify the present review. First, the corre- 
lators used as input are calculated now with greater preci- 
sion in perturbative QCD, in many cases up to the order 
af . Also, calculations with greater precision in Chiral Per- 
turbation Theory (ChPT), Heavy Quark Effective Theory 
(HQET), Soft CoUincar Effective Theory (SCET), or on 
the lattice, provide improved values of the form factors at 
some specific points. The techniques presented here are 
the optimal frame of including inputs coming from scpa- 
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rate sources and testing their consistency. Moreover, im- 
proved information about the phase of the form factor is 
available by Watson's theorem from the associated 
elastic scattering, and in some cases also the modulus is 
measured independently along a part of the cut. 

We first present in section [5] our notation and formal- 
ism and in section |3] the general Meiman interpolation 
problem, for an arbitrary number of derivatives at the 
origin and an arbitrary number of values at points inside 
the analyticity domain. The solution is obtained either 
by Lagrange multipliers, or by the techniques of analytic 
interpolation theory [571138] . and is written in two equiv- 
alent ways, as a determinant of a suitable matrix and as 
a compact convex quadratic form. 

In section 2] we present the complete treatment of the 
inclusion of the phase on a part of the unitarity cut, along 
with an arbitrary number of constraints of the Meiman 
type. No such treatment is found in the literature despite 
the long history of the problem [5.,6.9.26.,3IJ|. We derive 
two equivalent sets of integral equations, using, as in sec- 
tion [3J either Lagrange multipliers or the analytic inter- 
polation theory. 

In section [5] we treat the situation when, in addition to 
the phase, some information on the modulus along a part 
of the cut is available from an independent source. This 
is a mathematically more complicated problem [8lll2j. In 
this paper we shall present an approach proposed in [55], 
which uses the fact that the knowledge of the phase al- 
lows one to describe exactly the elastic cut of the form 
factor by means of the Omnes function. The problem is 
thus reduced to a standard Meiman problem on a larger 
analyticity domain. The method is very powerful and was 
recently employed to provide stringent constraints on the 
scalar Kit form factor at low energies |33j . However, while 
the method of treating the phase in section 2] automati- 
cally takes into account all the constraints of the original 
problem, this is not the case with the technique of this 
section: it provides necessary constraints for the input, 
which however may violate the original unitarity inequal- 
ity. Therefore, the allowed domain for the input values is 
given by the intersection of the domains derived by the 
formulas of sections H] and |S] We illustrate this fact in 
section |6l where we present constraints on the slope and 
curvature of the K13 scalar form factor. This section ex- 
tends the results previously reported in [23] • Finally, in 
the last section we summarize our conclusions and discuss 
possible applications. 

This paper provides a comprehensive treatment of im- 
portant mathematical and theoretical tools that are essen- 
tial for improving precision studies of form factors that are 
of great importance to the standard model. 



2 Notation and Formalism 

Let F(t) denote a generic form factor, which is real an- 
alytic in the complex t-plane cut along the positive real 
axis from the lowest unitarity branch point <+ to 00. The 
essential condition exploited in the present context is an 



inequality of the type: 



dt pit)\F{t)\-' < I, 



(1) 



where p{t) > is a positive semi-definite weight function 
and / is a known quantity. As mentioned in the Intro- 
duction, such inequalities can be obtained starting from a 
dispersion relation satisfied by a suitable correlator, eval- 
uated in the deep Euclidean region by perturbative QCD, 
and whose spectral function is bounded from below by a 
term involving the modulus squared of the relevant form 
factor. An example will be presented in section [S] 

In the analysis of the semileptonic decays one is in- 
terested in the parameters of the Taylor expansion at the 
origin, written as 



m = m 



2M4 



(2) 



where M is a suitable mass and A' and A" denote the 
dimcnsionless slope and curvature, respectively. Also, from 
the symmetries of QCD at low energies and, more recently, 
from lattice calculations, one may know F{t) at several 
special points inside the analyticity domain. The standard 
unitarity bounds exploit analyticity of the form factor and 
the inequality ([T]) in order to correlate in an optimal way 
these values and the expansion parameters in 

Additional information on the unitarity cut can be in- 
cluded in the formalism. According to Watson's theorem 
[5S] . below the inelastic threshold tin the phase of F{t) 
is equal (modulo tt) to the phase 6{t) of the associated 
elastic scattering process. Thus, 



F{t + ie) = \F{t)\e'^<^'\ 



t+ <t < U 



(3) 



where d{t) is known. Moreover, in certain cases also some 
information on the modulus or a bound on it, is 

available on the same range < t < tin ■ In section [3] 
we consider the standard version of the unitarity bounds, 
with no information about the phase and the modulus, 
except the inequality (|T]). The inclusion of the phase and 
modulus will be discussed in sections S] and [5j 

For the subsequent treatment, the problem is brought 
to a canonical form by making the conformal transforma- 
tion 

(4) 



that maps the cut i-plane onto the unit disc | z | < 1 in the 
z = z{t) plane, such that is mapped onto z = 1, the 
point at infinity to z = — 1 and the origin to z = 0. After 
this mapping, the inequality ([T]) is written as 



1- r\l%(e^^)p< J, 

Jo 



where the analytic function g{z) is defined as 
giz)^Fitiz))wiz). 



(5) 



(6) 
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Here t{z) is the inverse of (U) and w{z) is an outer function, 
i.e. a function analytic and without zeros in |z| < 1, such 
that its modulus on the boundary is related to p(0) and the 
Jacobian of the transformation (j?]). In general, an outer 
function is obtained from its modulus on the boundary by 
the integral 



w{z) = exp 



1 



2tt 



(7) 



In particular cases of physical interest, w{z) has a simple 
analytic form (for an example see section IH]). 

The function g{z) is analytic within the unit disc and 
can be expanded as: 



g{z) = go + giz + g2z^ + 
and (HI) implies 

CJO 



(8) 



(9) 



Using the real numbers gk are expressed in a straight- 
forward way in terms of the coefficients of the Taylor ex- 
pansion ([2]). The inequality represents the simplest 
"unitarity bound" for the shape parameters defined in ([5]). 
In what follows we shall improve it by including additional 
information on the form factor. 



3 Meiman problem 

We consider the general case when the first K derivatives 
of g{z) at z = and the values at N interior points are 
assumed to be known: 



1 d^-g(z) 



fc! dz*-' 



=0 

gizn) 



9k, 



0<k<K- 
l<n< N, 



1; 



(10) 



where gk and ^„ are given numbers. They are related, by 
jH), to the derivatives F^^\0), j < of F{t) at t = 0, 
and the values F{t{zn)), respectively. For simplicity and 
in view of phenomenological inputs that we will use, we 
assume the points z„ to be real, so ^„ are also real. 

Meiman problem |13| requires us to find the optimal 
constraints satisfied by the numbers defined in ((TO)) if 
(O holds. This mathematical problem is also known as 
a general Sehur-Caratheodory-Pick-Nevanlinna interpola- 
tion [Tl[15l[16]. 

One can prove [571155] that the most general constraint 
satisfied by the input values appearing in ()10[) is given by 
the inequality: 

Mo < I, (11) 
where /Iq is the solution of the minimization problem: 

(12) 



Here ||5||^2 denotes the norm, i.e. the quantity appear- 
ing in the l.h.s. of ([5]) or ([9|), and the minimum is taken 
over the class Q of analytic functions which satisfy the 
conditions (flTI)) . In the next subsections we shall solve the 
minimization problem (|12p by two different methods. 



3.1 Lagrange multipliers 

One may set up a Lagrangian for the minimization prob- 
lem (fT2l) with the constraints (fTOl): 



k=0 



N 

E 



(13) 



where a„ are Lagrange multipliers, and ^„ are known 
numbers defined as 



(14) 



fc=0 



Solving the Lagrange equations obtained by varying with 
respect to gk for all fc > and eliminating the Lagrange 
multipliers yields the solution of the minimization prob- 
lem (fT^ . For purposes of illustration, when N ~ 2, the 
Lagrange equations yield 



gk=aiz'[+a2Z^, k>K, (15) 
and the solution of the minimization (|12p takes the form: 

K-l 



"2 4 7 



Mo = E 5fe 

k=0 



CJO 

E 

k=K 



a2j29k4- (16) 

k=K 



Then inequality (jlip can be expressed in terms of the two 
Lagrange multipliers as: 



where 



aiCi + "26 < I, 

K~l 
k=0 

and the constraint conditions themselves are 



(17) 
(18) 



.,2K 



ai 



ai- 



l-zl 
(^1^2)^ 



Oi2 



{Z1Z2 



1 - ^12:2 



6, 



,2K 



Z1Z2 



OL2- 



(19) 



The consistency of eqs. ([T7)) and ([T^ can be written as: 



I2 



•2K 



6 

(^1^2)^ 



'% 1 
{Z1Z2) ^ 
1 - Z1Z2 1 



^1^2 
2K 
2 



> 0. 



(20) 



This can be readily extended to the case of N constraints: 



6 

^2K 
Zl 



6 

(^1^2)^ 



(ziZn)'' 



1 



z 



1 - Z1Z2 

2K 



6 



{ziZ2r {Z2) 

1 - ZIZ2 



(ziZn) {z2Zn) 



K 



1 - zizjv 1 - Z2ZN 



1 - ZlZN 
{Z2ZN)'' 

1 - Z2ZN 



^2K 

Zn 

1-4 



> 0. 



(21) 
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Alternatively, the solution can be obtained by introducing 
Lagrange multipliers also for the given coefficients gk, k = 
0, if — 1 in as was done in ref. [TU]. This leads to 
the inequality, equivalent to (j2ip : 



I 90 91 ■■■ 9K-1 

50 1 

51 f 



6 
1 

zi 



6 
f 

Z2 



1 

ZN 



I 

Zl 

2 



,K-1 



z. 



I — Z1Z2 1 — z 



I — Z1Z2 
1 



1-ZiZM 
1 



'TV 



1 — ZiZn 1 — Z2Zj^ 



> 0. 



OK-i 
6 1 Zl 
6 1 Z2 

1 ZN 

(22) 

This condition is expressed in a straightforward way in 
terms of the values of the form factor F{t) at ti ~ t{zi) 
and the derivatives at t = 0, using eqs. (0]) and (|6]). It can 
be shown that defines a convex domain in the space 
of the input parameters. 



3.2 Analytic interpolation theory 

Instead of using Lagrange multipliers, one can implement 
the constraints (jlOp by expanding the function 5 G ^7 in 
the most general way as [571(55] 



K-l 



N 



g{z) = QkZ^ + H ^^nBniz) + z^ B N +i{z)h{z) . 



k=0 



n=l 



(23) 

Here the functions Bn{z) are products of Blaschke factorfQ 
defined recurrently as 



Bi{z) = l,B„(z) 



Z - Zn-l 
1 - ZZn-l 



B„_i(z),2 < n<N+\, 



and the numbers An are obtained by solving the system 
of equations 



rri ^ 



K-\ 

E 

fe=0 



5fcZ„ 



1 < m < iV, 



where we took into account the fact that _B„(z,„) = 
for n > m. We recall that an expansion equivalent, but 
slightly different from ([23)). was proposed in [J. 

The function h[z) is analytic in \z\ < 1 and is free of 
constraints. Expressed in terms of ft-, the minimum norm 
problem ([T2]) becomes: 



{h} 



\h-H\\l2, 



(26) 



'^A Blaschke factor is a function B{z) analytic in \z\ < 1, 
which satisfies the conditions B{zo) = for a value \zo\ < 1 
and \B{e*'^)\ = 1 [14]. 



where iJ is a function defined on the boundary of the unit 
disc C = exp(i0) as 



1 



C^i3jv+i(C) 



K-l 



N 



LA:=0 



(27) 



and is meromorphic in \z\ < 1. 

The solution of ([26]) is straightforward. We expand: 



MC) = E'^'C', H{0^ E ^'C', (28) 



/ — — 00 



where Hi arc known real numbers defined as 



H, = 



1 



r'i/(c)^, 



-00 < / < 00. (29) 



ICI=i 



Using the expansions ([^5]) . we write (|26p as: 



Mo = min 



(30) 



where the minimization is taken upon the numbers hi. The 
minimum is reached for 



hi =Hi, l> 0, 

which leads to the minimum 

-1 



Mo 



= E^- 

/— — 00 



(31) 



(32) 



The coefficients Hi for I < —1 are calculated by inserting 
into ([29]) the expression of H from ([27]) and applying the 
residues theorem. The poles are produced by the factors 
Bn+i and for K + l-k> 0. The contribution of 

Bn+1 to iJ; is written as. 



TV 



(24) -J2 



BN+iiz) 



z—z. 



l-l 



K-l 



fc=i 



m— 1 



(33) 

where in the parenthesis we recognize from (|25p the num- 
bers S,nl z^ ■ The factor Q^+i-^ contributes to Hi as 



(25) -YO{K + l-k) — 



9k 



^^K+l-k 



k=0 



{K + 1 ^ ky. dz^+^-i' 



1 



Bn+i{z) 



z=0 

(34) 

Collecting the terms ([55]) and ([5^ we obtain, for / < — 1, 



N 



Hi = -Y^TkT-Y.^(^ + ^-^^ (35) 



K-l 



n=l " 

where we denoted 

Z Zri 



Yn 



hi = 



-B7V+l(-2) 
1 



fc=0 



z—Zj^ 

^K+l-k 



{K + I - k)\ dz^i+^-k 



1 



B]s[+i{z) 



(36) 



:=0 
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Using ([55)) it is easy to calculate the sum required in 
Due to the 9 function, only the values ^ > A: — in the 
second term of Hi give non-vanishing contributions. The 
result is written in a compact form as: 

TV K-l N K-1 

+ 2 ^ ^ Ckn9kS,n, 
j,k=0 n=l k=0 

(37) 



m,?i— 1 

where we defined 

Y Y 1 



" l=k^K 



E 



(38) 



and L = max(fc — K, j — K). Inserting (|37|) into the in- 
equality (|TT|) gives the allowed domain of the input values 
appearing in (fTU)) : 

N A'-l AT K-1 

y^^Amn£,n£.m + ^ BjkQjgk + 2^ ^ Ckn9k(,n < I- 
m.n—1 j^k—0 n—1 k—0 

(39) 

It can be checked that the domains given by ([?!]) and p9p 
are equivalent. 



4 Inclusion of the phase 

In this section we shall impose the condition that the 
phase of the form factor is known (modulo tt) along the 
elastic part of the cut from the phase of the associated 
elastic amplitude by Watson's theorem [36]. We start by 
defining the Omnes function 



dt 



(40) 



where S(t) is known for i < tin, and is an arbitrary func- 
tion, sufficiently smooth (i.e. Lipschitz continuous) for 
t > tin. From © and (gOl) it follows that 



Im 



F{t + ie) 
0{t + ie) 



0, 



t+ < t < tin. 



(41) 



Expressed in terms of the function g{z) this condition be- 
comes 



Ini 



WiO) 



0, 0e{-0in,Oin). 



(42) 



Here Sin is defined by z(tin) = exp(i0in) and the function 
W{9) is defined as: 



W{e) = w(e'^)0{e''>), 
where 10(2) is the outer function and 
0(z)=0(t(z)). 



(43) 
(44) 



As shown in [51 1^150] , the constraint (|^^ can be imposed 
by means of a generalized Lagrange multiplier. The con- 
straints at interior points can be treated either with La- 
grange multipliers as in subsection 13. 11 or by their explicit 
implementation as in subsection l3.2l Below we shall briefly 
present these two approaches. 



4.1 Lagrange multipliers 

The Lagrangian of the minimization problem (jl2p with 
the constraints (ITOl) and (1421) reads 



^ = ^ E -9' + E E akz") (45) 

k=Q n=l k=0 

00 ^i" 

+ -y5felim / Xi0')\Wi0')\lm[[W{e')]-\''e'''^]d9'. 



k=0 



The Lagrange multiplier \{0) is an odd function, A(— 6*) = 
—X{0) and, as in p51l30] . the factor |M^(6')| was introduced 
for convenience. 

We minimize C by brute force method with respect to 
the free parameters gk with k > K. The Lagrange mul- 
tipliers X{0) and a„ are found in the standard way by 
imposing the constraints pUj) and (|^^ . This leads to a 
system of coupled equations, which can be solved numer- 
ically. 

The calculations are straightforward (see for instance 
pS]). In order to write the equations in a simple form, it is 
convenient to define the phase 'P{0) of the function W{6) 

by 

Wi9) = \Wi0)\e"^'-''\ (46) 
From we have 



(47) 



where (f>{0) is the phase of the outer function w(e'*) and 
<5(t) is the clastic scattering phase shift. We introduce also 
the functions /3„ for n = 1, ...N, by 

o K Sm[K0-m]-MiK-l)0-m] 

= l + .2_2z„cos0 ■ 

Then the equations for the Lagrange multipliers X{0) and 
an take the form: 



K-l 



N 



k=0 



J2 9k sin[k0 - ^0)] = X{0) - J2 anPn{e) (49) 

d0'x{e'))ic^{e,e'), e e {-e,^,e,^), 



1 

'2^ 



"in N 

i j X{0)l3n{0)Ae+Y,^n\ 

a. 



{^Zji Zfi' ) 



K 



■^n, (50) 
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where n = 1, . . . N. 

The integral kernel in (|49l) . defined as 

_ sin[(j^ - l/2)(g - e') - m + 
^*(^'^) = sin[(0-0')/2] ' ^''^ 

is of Fredholm type if the phase ^{0) is Lipschitz continu- 
ous 130] • Then the above system can be solved numerically 
in a straightforward manner. Finally, the inequality (|11D 
takes the form: 



-. K-l N 

-Y.9k / d0A((?) sin [ke - ^{9)] + ^ a„f"„ < /, (52) 

fc=0 a "=1 



with / defined in ([18)) . Using the relation (|6]), the above 
inequality defines an allowed domain for the values of the 
form factor and its derivatives at the origin. Note that 
removing the phase constraint gives back the results of 
subsection 13.11 The results when phase alone is included 
{N = 0) as in [55] (and references therein) and the case 

= 1 discussed in 15011^ are also readily reproduced. It 
must be emphasized that the theory for arbitrary number 
of constraints is being presented here for the first time. 

It is easy to see that, if iin is increased, the allowed 
domain defined by the inequality (j52[) becomes smaller. 
The reason is that by increasing ijn the class of functions 
entering the minimization becomes gradually smaller, 
leading to a larger value for minimum /i§ entering the 
definition ([TT|) of the allowed domain. 



The second sum in the r.h.s. was already evaluated in 
subsection 13.21 The first term involves the function A(0), 
which we determine by imposing the condition (|42[) . writ- 
ten in terms of the function h{z) expanded as in (|55|) . 
Using hi given by ([55]) as the sum Hi + ci, we obtain by 
a straightforward calculation the integral equation for a 
function X(9) in the interval 9 £ {—9in,6in)'- 

fin 

A(6')-^ / d6''A(6'')/C^(6',6i') = ^(6*), (57) 

-fin 

where the kernel JC^ is defined as in ([51]) in terms of the 
known function 

'F{9) = arg[i3jv+i(cxp(i0))] - arg[T^(0)], (58) 

and \^ is a known function defined as: 

N 



Vie) = y: 

n=l 
K-l 

+ E 



r„Cn zn sin[if (0)] - sin[if (0) - 9] 



1 + — 2Zn COS ( 



9k 



k=0 



{K~k-iy. 



Uiz) 



(59) 



with 

Uiz) 



; sm[K9 + !f (6*)] - sin[(i^ - 1)61 + if (6*)] 
Sjv+i(z)(l + z2 -2zcos6') ■ 



(60) 



4.2 Analytic interpolation theory 

Alternatively, we shall implement the constraints in the 
interior points by expressing the function giz) as in eq. 
([55]) of subsection 13. 2[ in terms of a function /i(z) free of 
constraints. The Lagrangian can be expressed in terms of 
the coefficients hi defined in ([551) as: 



Using the expression ([M[) of ci and the integral equation 
([57]) . it is straightforward to evaluate the sum in the first 
term of ([5^ 



fin 

OO „ 

^c?^- / dOX{9)V{9). 

1=0 i. 



(61) 



C ^Y^^hi - Hi)'' + J2 Hf -2Y,hici + --- (53) 

1=0 



1=0 



/ — — OO 



where 



ci 



d6''A(6i')|M/'(6'')|Im 



eiO'(K+i)Bj,+,{e''') 
WiO') 



(54) 

The minimization of the Lagrangian given in ([531) with 
respect to hi has the simple solution 



hi = Hi+ ci, 
leading to the minimum 



> 



Mo 



(55) 



(56) 



Collecting all the terms in ([55|. ([TT]) can be written as: 



N K-l 
- / A9\i9)Vi9) + ^raninU + ^ Bjk 



9] 9k 



m.n— 1 
N K~l 



j,k=0 



1=0 



/— — OO 



+ 2 E E ^kn9k^n < I- (62) 
n=l k=0 

This inequality gives the allowed domain for the input 
values appearing in the conditions (|10p when the phase is 
known in the elastic region. We note that the arbitrary 
function Sit) for t > tin entering the Omncs function ([^Of 
does not appear in the result. The first term in ([55]) repre- 
sents the improvement brought by the information on the 
phase, as can be seen by comparing with ([39p . It can be 
checked numerically that the allowed domains described 
by ([52]) and ([62[) are equivalent. 
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5 Inclusion of phase and modulus 

In some cases, information on the modulus of the form 
factor along an interval of the unitarity cut is available 
from an independent source. As we mentioned in the In- 
troduction, a rigorous implementation of this information 
is difficult. In this section, we shall present an approach 
proposed in |26j , which leads to an independent constraint 
that should be satisfied by the inputs. 

We start with the remark that the knowledge of the 
phase was implemented in the previous section by the 
relation (HIJ, which says that the function f{t) defined 
through 

F{t) = /(t)O(t), (63) 

is real in the elastic region. In fact, since the Omnes func- 
tion 0{t) fully accounts for the elastic cut of the form 
factor, the function f(t) has a larger analyticity domain, 
namely the complex t-plane cut only for t > ti^. Imple- 
menting this fact leads to a modified version of the unitar- 
ity bounds, proposed in [55]. The method requires also the 
modulus of the form factor in the elastic region. Indeed, 
II]) implies that / satisfies the condition 



dtp{t)\0{t)\'\fit)\'<l' 



where 



dtpit)\Fit)\' 



(64) 



(65) 



If the modulus \F{t)\ is known for i+ < i < tin, the quan- 
tity /' is known. Then leads, through the techniques 
presented in section [31 to constraints on the values of / 
inside the analyticity domain. 

The problem is brought into the canonical form by the 
transformation 



(66) 



which maps the complex t-plane cut for t > tin on to the 
unit disc in the z-plane defined by z — z{t). Then ([64)) can 
be written as 



^£^d%(cxp(z0))P</', 



where the function g is now 

g{z) = w{z)u:{z)F{i{z))[0{z)]-\ 



(67) 



(68) 



Here w^z) is the outer function related to the weight p{t) 
and the Jacobian of the new mapping ([55)) and 0{z) is 
defined as 

0{z) = 0{i{z)), (69) 

where t{z) is the inverse of z = z{t) with z{t) defined in 
([55)) . and 



uj{z) = exp 



V^tin - t(z) 



dt'— 



\n\0{t')\ 



Unit' ~ t(z)) 



The inequality ([57)) has exactly the same form as ([5]). 
Therefore, by using the techniques of section [3| we derive 
constraints on the function g at interior points, which can 
be written in the equivalent forms as in ([?!)) or ()39p . Us- 
ing ([55)) . these constaints are expressed in terms of the 
physically interesting values of the form factor F{t). 

In fact, the Omnes function 0(t) defined in ([40|) is not 
unique, as it involves the arbitrary function 5{t) for t > tin- 
We have seen that the results of section [J] are not affected 
by this arbitrariness. This is true also for the results of 
this section: the reason is that a change of the function 
S{t) for t > tin is equivalent with a multiplication of g{z) 
by a function analytic and without zeros in \z\ < 1 (i.e. 
an outer function). According to the general theory of an- 
alytic functions of Hardy class [H1[T5] , the multiplication 
by an outer function does not change the class of functions 
used in minimization problems. In our case, the arbitrary 
function 6{t) for t > tin enters in both the functions 0{z) 
and uj{z) appearing in (|68p. and their ambiguities compen- 
sate each other exactly. The independence of the results 
on the choice of the phase for t > tin is confirmed numer- 
ically, for functions 6{t) that are Lipschitz continuous. It 
is important to emphasize that the method relies on the 
Omnes function making an appearance first through a re- 
lated outer function and then through its inverse, while 
the function /(t) is merely introduced at an intermediate 
stage and is subsequently eliminated. 

The constraints provided by the technique of this sec- 
tion are expected to be quite strong since they result from 
a minimization on a restricted class of analytic functions, 
where the second Riemann sheet of the form factor is de- 
scribed explicitly by the Omnes function. On the other 
hand, it is easy to see that the fulfillment of the condition 
([5^ does not automatically imply that the original condi- 
tion ([ij is satisfied. More exactly, the technique described 
here does not impose the knowledge of the modulus in 
addition and simultaneously with the bound ([T|) and the 
knowledge of the phase, but exploits only a consequence 
of the original conditions of the problem. Therefore, one 
must calculate separately the allowed domains of the pa- 
rameters of interest given by the techniques of sections [J] 
and[Sl and take as the final results the intersection of these 
domains. 



6 Example: scalar Kn form factor 

We consider as an example the scalar Ktt form factor 
/o(t), presenting constraints on the slope Aq and curva- 
ture Aq, appearing in the expansion 



./o(t) = /o(0) 



t 



t^ 



(71) 



(70) 



often used in the physical range of the semileptonic decay 
K — > ttIv. We work in the isospin limits, adopting the 
convention that Mk and Af^r are the masses of the charged 
mesons. 

The scalar Ktt form factor has been calculated at low 
energies in ChPT and on the lattice (for recent reviews see 
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[g^HOl HT]). At t = 0, the current value /o(0) = 0.962 ± 
0.004 [4^ shows that the corrections to the Ademollo- 
Gatto theorem are quite smalL Other low energy theorems 
frequently used are 



2.0 



MA 



Fk 



Fk 



-A 



CT, 



(72) 



where Ak-k = M]^ — and Axn = —Ak-k are the first 
and second Callan-Treiman points, denoted below as CTi 
and CT2, respectively. The lowest order values are known 
fr-om Fk/F^ = 1.193 ± 0.006 gD], and the corrections 
calculated to one loop are Act = —3.1 x 10"'^ and Act = 
0.03 [13|. In the isospin limit, the higher corrections are 
negligible at the first point, but are expected to be quite 
large at the second one. 

An inequality of the type ([1]) is obtained for /o(i) start- 
ing with a dispersion relation satisfied by a suitable cor- 
relator of the strangeness-changing current [30ll3T] : 



(73) 



where unitarity implies the inequality: 



Im7To(0 > - 



3M- [{t-t+){t-t^)]^/'- 



2 167r 



<3 



\foit)\\ (74) 



with t± = {Mk ± Af^)^. The quantity Xo(Q^) was calcu- 
lated up to the order in perturbative QCD [451. From 
([75]) and ([7^ it follows that in this case the quantity / 
appearing in ([5]) is 



xo(Q'). 



(75) 



For illustration we give also the outer function w{z) en- 
tering (O: 



w{z) 



32^^ Mk 



(l-z)(l + z)3/2 



f)f (l-Zz(t_))V2 



(l-zz{-Q'^)Y (H-z(i_))i/2 



(76) 



Further, below the inelastic threshold the phase is known 
from the J = 1/2 5'- wave of Kn elastic scattering 
while the modulus |/o(i)| was measured recently from the 
decay r — >■ KttVt- [57] . For details of the input quantities 
see [33] ■ When the modulus is also included, as in section 
([5]), the outer functions in ((68)) can be written as. 



Vi{Ml ~ Ml) ./T—z (1 + zf/\l + ~z{-Q^)f 



16y2^ti„ (l-zi(-Q2))2 
{I- z~z{t+)Y'^{l- z~z(t-)Y'^ 

(l + z(i+))l/2(l + ;5(t_))l/2 • 



(77) 



In ref. |33| we derived stringent bounds on the slope 
Aq and curvature Ag, using as input the values of fo{t) at 
t = and CTi, and information on the phase and modu- 
lus included in the formalism of section [S] The additional 
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Fig. 1. Allowed regions in the slope-curvature plane of the 
scalar Kit form factor, obtained with the techniques of sections 
00 For explanations see the text. 



input from the unitarity cut led to a dramatic improve- 
ment of the bounds obtained in [53] with the standard 
unitarity bounds of section |31 In [33] we obtained also a 
narrow allowed range for the ChPT correction Act at the 
second Callan-Treiman point CT2. 

In the present paper we further illustrate the tech- 
niques presented in sections |3l |4] and |5j by comparing 
their constraining power for various inputs. In fig. |T] we 
present the allowed domain for the slope and curvature 
of /o(i), obtained with the methods described above: the 
large ellipse is obtained with standard unitarity bounds of 
section[31 the intermediate ellipse includes the phase up to 
<in = 1.0 GeV^ with the method of subsection 14. 1[ and the 
small ellipse is obtained with the method of section jS] for 
the same tin. The domains shown in the left panel are de- 
rived using as input the normalization /o(0) = 0.962 and 
/o(A;f7r) = 1.193, those in the right panel use as input also 
a value at the second Callan-Treiman point f^{AK-n) = 
1/1.193 + Act for a certain choice Act = -6.0134. 

As emphasized in section |4l the domain obtained with 
the inclusion of the phase should be contained entirely 
inside the domain obtained with the expressions of section 
131 since it is related to a minimization on a smaller class 
of admissible functions. This is confirmed by the large and 
intermediate ellipses of fig.|T] On the other hand, the small 
ellipses are given by a minimization on an admissible class, 
defined in section [SJ which is not a priori contained in the 
class defined in section|4l Therefore, the small ellipses need 
not to be contained entirely in the intermediate ones. 

In the left panel of fig. |TJ the small ellipse is contained 
inside the other two, which means that for points inside 
this allowed domain all the constraints are fulfilled. In the 
right panel, where we use as input also the value at CT2 , 
all the domains shrink, and the small ellipse has a part 
situated outside the intermediate one. To satisfy all the 
constraints, one should take the intersection of the small 
and intermediate ellipses. 
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Clearly, by increasing ijn, i.e. the energy up to which 
the phase (and the modulus) are given, the system is more 
and more constrained and one may reach a situation when 
the inputs become inconsistent. This is illustrated in figs. 

and [31 where we show the configuration of the allowed 
domains for two larger values of ti^. The large ellipses are 
the same in all figures, since they are are independent of 
tin- As follows from the arguments of section SI the in- 
termediate ellipses become gradually smaller when tin is 
increased. In the left panels, obtained with the normal- 
ization at t = and the value at CTi, the small ellipses 
are contained inside the intermediate ones, indicating that 
we can find an allowed domain that satisfies all the con- 
straints. However, if we impose also the constraint at CT2, 
the system becomes over-constrained, and it is impossible 
to find a domain that satisfies all the constraints. Indeed, 
in the right panels of figs. [5] and |31 the small ellipses are 
not inside the intermediate ones. 
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Fig. 2. As in fig.[T]for tin = (1.2 GeV)^. The inset zooms in on 
the smallest ellipse obtained using the techniques of section [S] 



The same behavior is illustrated in figs. 2] and El where 
we show the allowed regions in the slope-curvature plane 
for various values of the inelastic threshold t\^. The left 
panels show that by increasing tin the ellipses obtained 
with the phase included as in section 21 become gradually 
smaller and are contained one within the other. This is 
not the case with the ellipses derived with the phase and 
modulus included as in section |S1 In the particular case 
when both CTi and CT2 are used as input, the ellipses 
obtained with various tin have a zero intersection (see right 
panel of fig. [5]). This signals an inconsistency of the various 
quantities used as input (normalization at t = 0, values 
at CTi and CT2, the phase up to tin and the modulus 
below tin). In particular, the assumption of neglecting the 
inelastic effects up to tin = (1.4 GeV)^ may be too strong. 
We note however that the above figures were obtained 
with the central values of the inputs described in |33j . Of 
course, in phenomenological analyses one should account 
for the errors of the various pieces of the input. We did 
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Fig. 3. As in fig.[T]for tin = (1.4GeV)^. The inset zooms in on 
the smallest ellipse obtained using the techniques of section (5] 



not consider this aspect here, because our purpose was 
only to illustrate the mathematical techniques in a definite 
framework. Thus we have provided a concrete illustration 
of the various techniques discussed in detail in the previous 
sections. The phenomenological implications of the results 
will be analyzed in a future work. 
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Fig. 4. Allowed region in the slope-curvature plane for various 
tin using as input the normalization at t = and the value at 
CTi. Left: phase included as in section (4] Right: phase and 
modulus included as in section [S] 



7 Conclusions 

In this paper we reviewed the method of unitarity bounds, 
extended in order to include information on the phase and 
modulus of the form factor on the unitarity cut. We pro- 
vided explicit formulas, easily implementable in Mathe- 
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Fig. 5. Allowed region in the slope-curvature plane for various 
tin using as input the normalization at t = and the values 
at CTi and CT2. Left: phase included as in section [l] Right: 
phase and modulus included as in section [S] 



matica or C programs, for an arbitrary number of deriva- 
tives at t = and an arbitrary number of interior points. 

The method is very suitable for correlating through 
analytieity various pieces of information about the form 
factors: perturbative QCD, lattice calculations and effec- 
tive theories like ChPT or SCET. It does not depend on 
specific assumptions usually adopted in the standard dis- 
persion relations, like the absence of zeros, or the behavior 
of the form factor above the inelastic threshold. 

As shown in section |6l the techniques presented here 
provide strong constraints on the shape parameters of 
the form factors, of interest for the parameterizations of 
the experimental data. The method can be used also to 
test the low energy theorems and put constraints on the 
higher order corrections of ChPT [Sj. Moreover, it can 
be adapted in order to control theoretically the trunca- 
tion error of the experimental parametrization |251I35| . 

Another possible application is the detection of the ze- 
ros of the form factors. The problem is of interest for the 
Omnes representations used recently for the parametriza- 
tion of various form factors, which assume that the ze- 
ros are absent. For such a study, one has to assume that 
Pito) = for a certain unknown to, include this condition 
among the interior constraints ([TU)) and test the consis- 
tency of the inputs, in one of the versions presented in 
sections [3][5l The method then can give in an unambigu- 
ous way the points to where zeros are excluded. 

The method of unitarity bounds proved to be very use- 
ful for the description of the B^D,B^tt 01 K^tt 
form factors. It can be applied also to other form fac- 
tors, such as those describing D — >■ tt semileptonie decays, 
or the scalar form factors of the pion or kaon, for which 
bounds of the type H]) can be obtained from the dispersion 
relation of a suitable correlator calculated in perturbative 
QCD. The method is a valuable tool for increasing the 
precision of the predictions in low energy flavor physics. 



which has been discussed here in great detail and gener- 
ality in an accessible manner. 
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